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Abstract 

The question of how infalhng matter in a pure state forms 
a Schwarzschild black hole that appears to be at non-zero temperature 
is discussed in the context of the AdS/CFT connection. It is argued 
that the phenomenon of self-thermalization in non-linear (chaotic) sys- 
tems can be invoked to explain how the boundary theory, initially at 
zero temperature self thermalizes and acquires a finite temperature. 
Yang-Mills theory is known to be chaotic (classically) and the imagi- 
nary part of the gluon self-energy (damping rate of the gluon plasma) 
is expected to give the Lyapunov exponent. We explain how the imag- 
inary part would arise in the corresponding supergravity calculation 
due to absorption at the horizon of the black hole. 
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1 Introduction 



Ever since the discovery by Bekenstein and Hawking that black holes have 
an entropy [Q and a temperature [§] associated with them there has been a 
puzzling question. When infalling matter in a pure quantum state collapse 
into a black hole, how does it become a mixed state, which a Schwarzschild 
black hole with a finite Hawking temperature must needs be? One possible 
answer is that quantum mechanics in the presence of gravity is not unitary. 
Another is that the state is actually pure and only appears mixed to the 
outside observer. In this paper we wish to dwell not on this question but 
rather on a related simpler question: Why does the state appear thermal to 
the outside observer (regardless of whether it is actually pure or mixed)? This 
second question in some sense can be phrased in an operational manner since 
a thermal state (in the thermodynamic sense of the word) can be defined to 
be one which for instance has a thermal (Maxwell-Boltzmann, Bose-Einstein 
or Fermi-Dirac) distribution of particles. Or it can be defined as a system 
that looks ergodic. 

Chaotic systems typically become ergodic on a time scale given by the 
(inverse of the) Lyapunov exponent. Thus it is possible for a non-thermal 
state to become thermal. Therefore an answer to the above question could be 
that the system is chaotic. Thus we could say that infalling matter initially 
in a coherent non thermal state must self thermalize upon collapsing into a 
black hole. 

Of course this problem must be treated quantum mechanically. Quan- 
tum chaos is less well understood, but if we accept Berry's conjecture as a 
definition of quantum chaos then it has been shown that this thermalization 
takes place quantum mechanically also 0]. In fact these ideas can be used to 
define statistical mechanics without invoking a God-given "heat bath" 0. 

Of course it is rather difficult to study self-thermalization in a compli- 
cated situation such as black hole formation. The recent connection between 
gravity and boundary conformal field theory |^, |^, ^ allows us to pose 
this question in terms of fiat space Yang-Mills theory. It is thus reassuring 
to know that classical Yang-Mills exhibits chaos and self thermalization and 
that the Lyapunov exponent is nothing but the damping rate of hot (quan- 
tum) gluon plasma. This has been calculated for high temperature QCD 
[§, ^, |Ty] and seems to agree with numerical results on chaos in Yang-Mills 
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What is found [|I^ is that the damping rate of a thermal plasmon at rest 
is obtained in perturbation theory to be [Q: 

70 ^ 6.635— g'T 
= 0.176/r = 2 

= 0.26VT = 3 (1.1) 
The Lyapunov exponent on the other hand is found numerically to be 

MM 

Ao ^ OMg^T {N = 2) 

^0.53/T (A^ = 3) 

It is thus seen that ^= 2. (As explained in |]10| the factor of two can be 
traced to the way the two quantities are defined.) The numerical claculations 
[|ri[] are done on a lattice of approximately (20)^ sites and for a range of values 
of the lattice spacing and g. It is found to be chaotic even for very small 
values of g which measures the non-linearity. At very high temperatures 
one might expect that classical results have some validity in the quantum 
regime also. This is also reflected in Berry's conjecture about the nature 
of quantum energy levels of a classically chaotic system - the assertions are 
about the higher energy levels. Therefor we assume that QCD with large 
enough energy density exhibits thermalization. 

If we further make the (reasonable) assumption that at high temperatures 
Supersymmetric Yang-Mills is qualitatively similar to QCD (at high temper- 
ature) the puzzle is resolved: When we add some coherent energy density 
to a Yang Mills vacuum the system thermalizes and becomes ergodic and 
then one can associate a temperature T related to the energy density in the 
usual way E/V ~ T^. In gravity this process corresponds to the formation 
of a black hole with a finite Hawking temperature. This is an answer to 
the ("simpler") question posed above at least at a conceptual level. It still 
remains to actually calculate this exponent in supergravity. This is similar 
to the calculation of the glueball mass using supergravity. 

The real part of the self energy in the confining phase, the glueball mass, 
in 3+1-QCD was calculated using supergravity |T2|, |T3|, [I^. Similarly one 
can calculate the imaginary part of the self energy in high temperature 3-1-1- 
Supersymmetric Yang-Mills. Let us understand the supergravity origin of the 
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imaginary part. From simple arguments one can show that the eigenvalues 
of a Schroedinger like operator have imaginary parts when the boundary 
conditions are such that the current going into a region is not equal to the 
current going out. If we imagine say a (normalizable) solution describing a 
dilaton outside the horizon of a black hole, with incoming flux equalling zero 
at infinity, then, as the dilaton falls into the black hole the flux at the horizon 
will be non zero. Such a solution must have an imaginary part. It is crucial 
here that the boundary condition just outside the horizon is that there are 
only ingoing waves. These are complex boundary conditions and that is why 
the eigenvalues are complex. In Appendix A we give a very brief description 
of a one-dimensional quantum mechanical problem that illustrates this point. 
In the glueball mass calculation of [0, jT^], |T^ |T3[, the boundary condition was 
that the derivative vanish at the horizon. This ensures that the eigenvalues 
are real. This boundary condition is correct for an equilibrium situation. 
Another way to see that this is correct is to look at the imaginary time 
description of equlibrium statistical mechanics. In this case the horizon is 
like the origin of a polar coordinate system and consisitency demands such 
a boundary condition. The non-equlibrium situation on the other hand does 
not require such a boundary condition. It requires the complex boundary 
condition that we described above. 

One can then ask: If the complex boundary conditions corresponding 
to non-equilibrium are due to absorption by a black hole, how can there 
ever be equilibrium in such a situation i.e. how can a dilaton possibly not 
be absorbed by a black hole? The answer, obtainable only in a complete 
theory, is likely to be that Hawking radiation takes care of the emission and 
in equilibrium absorption and emission must be equal. The imaginary part 
is a measure of how fast equilibrium is reached after a small perturbation 
disturbs the system. It is equal to the difference between the two rates. In 
equilibrium the imaginary part will be zero. 

What does Hawking radiation correspond to in Yang-Mills? As it repre- 
sents a flow of energy from small to large radii, it must represent the flow of 
energy from long wavelengths to short wavelengths of the Yang-Mills fields, 
as the energy is equipartitioned in reaching thermal equilibrium. This is due 
to the same non-linearity of the equations that is also responsible for the 
phenomenon of chaos and self thermalization. In thermal equlibrium energy 
is uniformly distributed in all modes as demanded by equipartition. 

If the above idea is basically correct one can ask further questions, such 
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as what is the description within Yang-Mills of the phenomenon of a particles 
falling into a black hole. The answer is known A particle/wave packet 
in the bulk starting a distance 6 > away from the boundary corresponds in 
the boundary to a wave packet of width AX ~ 6. As the particle in the bulk 
falls radially inwards, AX increases, by the usual IR/UV duality between the 
bulk and the boundary. When the particle reaches the horizon, AX reaches 
the thermal wavelength. 

Also, the speed of the infalling particle cannot exceed that of light. This 
translates into an upperbound on the rate of wavepacket spreading in the 
boundary. This is causality in the bulk, and its manifestations in the bound- 



ary theory have been studied in ||T6 



One can also ask more refined questions such as what happens if the 
impact parameter of the bulk particle is sufficiently large so that the particle 
never reaches the horizon. When the impact parameter is non zero, the 
particle has angular momentum which is conserved in the bulk, certainly 
so in the supergravity approximation used to describe the bulk space time. 
Thus, as the particle falls inwards, its radial velocity decreases, its angular 
velocity increases. If the angular momentum is sufficiently large, then the 
particle will not reach the horizon at all. This would correspond in the 
boundary theory to a wave packet that never thermalizes. This is somewhat 
counterintuitive. 

Furthermore, the bulk angular velocity corresponds to the boundary wave 
packet momentum pi, in the boundary space time. The bulk angular mo- 
mentum conservation should then translate into the boundary theory as the 
conservation of the ratio of the wave packet momentum and its size. More- 
over, as the size increases (bulk particle's radial location decreases), the wave 
packet momentum increases ! It is not obvious that the boundary theory 
(Yang-Mills theory at finite temeprature) has such an unusual conserved 
quantity. 

It is in fact quite possible that there is no such conserved quantity. 
AdS/CFT corresondence would then predict that the angular momentum 
of the bulk particle is not conserved. This is not such an outrageous predic- 
tion. Presumably what is happening is that the bulk particle is interacting 
with the Hawking particles emitted by the black holes, thereby loosing its 
angular momentum. In this situation any particle will fall into the black 
hole eventually regardless of its initial angular momentum. While the super- 
gravity theory used for the bulk does not know of this effect, the boundary 
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theory somehow knows about this and is predicting its effect. One would 
expect that if the supergravity equations are corrected to include self inter- 
action of dilatons, equivalently the back reaction of Hawking photons, this 
effect would be seen. 

Regardless of whether one neglects the drag effect of the heat bath or 
not, one thing is clear from the above. If the wavepacket has a large value 
of the boundary momentum pi it will take longer to thermalize. It turns 
out that this is seen in numerical simulations for Yang-Mills - the Lyapunov 
exponents are found to decrease with pi. Thus the supergravity connection 
provides an intuitive understanding of this fact. 

If the ideas described above are correct one should be able to calculate the 
Lyapunov exponent for Yang-Mills using the supergravity connection. One 
should then compare with that obtained using finite temperature methods 
in QCD. One should of course keep in mind that the latter is computed 
in the perturbative regime whereas the supergravity calculation necessarily 
describes the strong coupling region. The inverse gluon two point function 
at finite temperature has the form (ignoring the index structure): 

where 11 (a;) contains the finite temperature corrections to the tree level result. 
This is calculated in the imaginary time formalism where u takes values 
27TinT. We can analytically continue u to the real axis and locate the poles 
^^0 + ^7o- Clearly 70 ~ ^"2^^^°^ • ^0 in fact is the damping coefficient. On the 
supergravity side rather than calculate the propagator, we determine uo + ijo 
directly by solving the dilaton equation. [] 

The computation proceeds as in except for two differences: One is that 
we will assume a real time dependence and solve for the energy eigenvalue 
in terms of the momentum. The second is that the boundary condition 
will be different and will lead to a complex eigenvalue. We will not be 
able to evaluate the eigenvalue analytically but we can argue that it will be 
complex and also determine the overall scale. The numerical evaluation will 
be reported elsewhere. 

^The real-time method does not reproduce correctly the imaginary part of the two point 
fmiction jl7[ |l8|. Nevertheless the location of the pole will be correct. This implies that 
the imaginary part of the self energy function which determines the Lyapunov exponent 
is also correct. 
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Finally we would like to make some comments on the issue of why classical 
Lyapunov exponent has any significance in quantum Yang-Mills. What we 
have calculated is the damping coefficient of a hot (quantum) gluon plasma 
- this is a classical calculation in supergravity, but it is fully quantum me- 
chanical in Yang-Mills. It turns out that in perturbation theory the leading 
term is independent of Ti (see eqn ( |1.1| )) and is a classical result |jlO|, ||. It 
also turns out, as mentioned earlier, that this leading term is equal to the 
Lyapunov exponent of classical Yang-Mills. Thus we conclude from this that 
classical Yang-Mills is chaotic. Now when we go to the quantum theory, 
presumably we must retain the h corrections to the gluon plasma damping 
coefficient. But the direct relevance of this quantity to quantum chaos is 
not known. As mentioned earlier, it is generally believed that systems that 
are classically chaotic are also quantum mechanically "chaotic" . If these are 
true then quantum Yang- Mills (with sufficicient energy density) will also look 
thermal. It is this property of classically chaotic systems that we are invoking 
to make statements about the quantum theory. 

To summarize, what we have calculated is the damping coefficient of 
hot (quantum mechanical) gluon plasma. This is a physically interesting 
quantity by itself. Moreover, its classical limit gives the Lyapunov exponent 
of classical Yang-Mills which, is thus chaotic. Qualitatively we expect that 
quantum Yang-Mills will also be chaotic and hence look thermal. 

The role of the h corrections (that we have in effect computed) to quan- 
tum chaos in quantum Yang-Mills is an open question. One can reasonably 
conjecture that if one had an exact effective action for hot (large N) QCD 
(i.e. after doing the functional integral) its Lyapunov exponent would be this 
fully quantum mechanical damping coefficient that we have computed. 



2 Preliminaries 

The metric describing an AdS^ - Schwarzschild black hole space-time is 
ds^ = - TIT + 1 :r)dt + —2 7777- + r^dn\ 

dr^ 

= -V{r)dt^ + j^ + r^dn^ (2.1) 
V[r) 
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The boundary is at r — oo. So we let the boundary be at r = i? and 
assume R is very large. The boundary metric is then 

= -dr^ + [dx'f (2.2) 

The period of t is 

In the above we have used r+ which is defined by: 

^W = ^(^'-^+)(^' + ^') 
Using — = — 6^, we get another expression for Ph 

Ph = T^— 4t 2.4 

Ph defines a Hawking temperature Th = The boundary temperature 
for the metric ( |2.2| ) is then given by Tb = ^Th- 

We can consider two limits. If we let GM << 6^ we get Tb = 2hi(^/^)- 

If we take the opposite limit we get Tb = (^P^) * • Note that the latter 
limit is used in along with a further change of coordinates ^ = (^p^)^. In 
this limit the boundary becomes as against for finite 



GM 

b2 ■ 



3 Classical Trajectories and WKB Solution 

If we let po = E, the energy, and L be the angular momentum, these being 
conserved, they serve to label the trajectories. The mass shell condition 
g'^'^Pij.Pu = gives us the equation 

1 ,dr,. L"^ , , 

+ T77t(3v) +^ = (3.5) 



V{r) V{r)'dX' r 
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Here A is an affine parameter that is defined by = Thus we get 

i^^r = E'-V{r)^ = E^-U{r). (3.6) 

The quantity U{r) is analogous to a classical potential and is shown in 
figure 1 in the case {^^) — > cxo. 




I 

! 

i Fig 1 , The classical "potential " when M » b, 
In rescaled coordmates where b is held finite, 



the hump m the potential is at infmity, 



The height of the potential is -^-(1 + 4^^?)- Thus whenever E'^ exceeds 
this the classical particle will fall into the black hole. Note also that in the 
limit taken in the hump is at infinity. So the region being dealt with is to 
the left of the hump. In this case all particles fall into the black hole. 

To get a WKB type solution we assume a solution of the form 

/ = e'*' (3.7) 

with 5* = 5*0 + ^5*1 + .... On substituting in the wave equation (see ( |4.9| ) 
below) one finds 

For large r this becomes: 

So = --^EW - L2 
r 
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3? 

Si = -lnr (3.8) 

The higher terms for large r can be seen to involve higher powers of ^^/gip^x^- 
Thus it is a good approximation for r « h and also far from the turning 
point. The turning point r* is when E"^ = U{r*). For r < r* the solution is 
damped, rather than oscillatory. Thus for r << r* the field vanishes. The 
flux thus vanishes and the dilaton is not absorbed by the black hole. As 
shown in section 5 one can show in such cases that the eigenvalue E is real 
and therefore the Lyapunov exponent is zero. Thus we see as claimed in the 
introduction, that in this approximation when the impact parameter is larger 
than the black hole radius the Lyapunov exponent is zero. 

4 Dilaton Wave Equation 

In the above metric the wave equation for the dilaton is: 

Here the dilaton field has been written as: 

$(t,r,(])=e-^V(r)lL(^^) 

where L is the angular momentum and Q represents all the angles, and s = lE 
(with E not necessarily real) specifies the time dependence. To be precise 
we can use the notation 

s = iuj + -f. (4.10) 

We will assume that uj,'-f are greater than zero. Note that 7 is the damping 
coefficient that we want to determine. 

We look for normalizable solutions. This equation ([4.9|) can also be writ- 
ten as 

^^2^ ^3v[r) dr V{ry V{r) r^^^ ^ ^ ' 

For large r as in the solutions either go to a constant which is not 
normalizable or go as ^ which is normalizable. Near r = r+ the equation 
becomes: 
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+ 4- + [7 ^ + 7 7]/ = 0- (4-12) 



dr'^ r — dr (r — r+)^ (r — r+) 

where A = — ^t^2 and 5 = —75 — f ^''^ ^ ^ 
Note that if we use (|2.4| ) we can rewrite 



A = (4.13) 



The substitution f = [r - r+)" E^I^T an(r - r+)" into ( ^l2|) gives the 



condition = —A. Thus we can refer to ±a as the two eigenvalues. We 
should also note that if we set the energy parameter s to zero we get degener- 
ate solutions. If we let y±a denote the two solutions, then in the degenerate 
limit a — i> 0, ^+°~^-" is the solution to consider and in fact goes as /n(r — r+). 

For a 7^ 0, a certain linear combination of these two y^a + P{c(, B)y+a 
is the one that is normalizable at infinity. For a given value of a and B 
there is a number (3{a,B) (that can be determined numerically) such that 
Va + (^{oi-i B)y_a is normalizable, i.e. the condition of normalizability fixes 
the functional dependence of /? Qon a and B. 

In the case at hand we will see below that the boundary condition at the 
horizon requires that /3(a, B) = 0. This is an eigenvalue condition and gives 
us the energy eigenvalues as a function of B or L. In the limit that — >■ 00, 

can be replaced upto some overall scale by /c^ the three momentum squared 



5 Boundary Conditions and Complexity of 
Eigenvalues 

We need to impose some boundary condition to fix the parameter f3 of the 
previous section. The boundary can be chosen to be at any value of r, 
and the boundary condition should reflect some physical requirement. As 
explained in the introduction, in [j^ the boundary condition was that at 
r = the derivative of the dilaton field vanish. The justification for this is 
that in equilibrium one should consider a Euclideanized time coordinate r as 

^Not to be confused with any inverse temperature. 
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is common in finite temperature field theory. Then the point r = r"*" is hke 
the origin of a polar coordinate system in the r — r plane. If the field is to 
be well defined the radial derivative must vanish. 

In our case we have a non-equilibrium situation. We have a real time 
coordinate. So this argument is no longer applicable. The physical problem 
we are trying to describe is that of a dilaton outside the black hole falling 
into the black hole. A reasonable boundary condition is that just outside 
the horizon the field must be purely ingoing. There cannot be any outgoing 
wave at the horizon. This is in accord with physical intuition that says 
once a particle reaches the horizon it cannot go out. We will thus impose 
this boundary condition at an infinitesimal distance /x outside the horizon at 
r = r+ + /i. /i is just a regulator that will be set to zero in the end. 

Let us first show that the eigenvalues have to be complex. We multiply 
by /* and integrate equation ( |4.!j| ) from r = a to r = 6 to get 



d 



Integrate by parts to get 



V{r) 



-]rfdr = 0. 



(5.14) 





df 


J a 


dr 



r^dr+[f*^r^V(r)] 
dr 



r=b 



+ 



^^^]rfdr = 0. (5.15) 
V{r) r^ 



Subtract from the LHS its complex conjugate. The first term is manifestly 
real and drops out of the difference. The second term minus its complex 
conjugate gives the difference in flux at the points r = a and r = b. In the 
last term only can possibly be complex and so we finally obtain 



[r^r^V{r)-c.c.] 



- (s 



* \ 2] J"* 



]rfdr 



(5.16) 



la V{r)' 

Let us take a = oo and 6 = r+ + yU. If we have a normalizable function / 
that vanishes at infinity and boundary condition at r+ + /i described above, 
the LHS is clearly non-zero. Thus s ^ s*. This proves that the eigenvalues 
are complex. 
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6 Solutions 



The two solutions to are 

n=oo 

y±a= E a„(r-r+r±" (6.17) 

n=0 

If one calculates the current ■^\/9g^^{4'*^ ~ c-c.) one finds that near 
r = r_|_ it is proportional to ±/m a. We will adopt the convention that 
Im a > 0. We want the current to be in the — r direction for r = r+. This 
means we pick y^a as the near horizon solution. Thus if the solution that is 
normalizable at oo is 

y{r) = y_a + Py+a : r > r+ + , (6.18) 

then the boundary condition requires that (3 = 0. 

This is an eigenvalue equation for a which we have seen has to have 
complex solutions. The actual value can be determined numerically. An im- 
portant point to note here is that this equation is independent of fi. Thus 
the Lyapunov exponent which is proportional to the real part of a is inde- 
pendent of fi. As described in this exponent is equal to the gluon plasma 
damping coefficient calculated in ^j. It was shown there that to leading 
order in k"^ (our B or L^) this quantity is a pure number times g^T in pertur- 
bation theory. In particular it does not have any infrared divergence. This 
last fact is presumably refiected in the /i-independence mentioned above. In 
fact we find that it is //-independent even for finite B. Thus the infrared 
divergences that were present in the perturbative QCD calculations Q are 
absent here. This is not surprising. In this formalism (which corresponds to 
strong coupling QCD) confinement and a mass gap are manifest. Therefore 
one does not expect infrared divergences. 

Also if we use ( |4.13| ) we see that a = Thus a being some pure num- 
ber independent of ratios of dimensionful parameters, we automatically find 
that s and hence 7 (see ( [4. 10] )) the damping coefficient, is some pure number 
times the temperature. We do not expect to reproduce the ^f-dependence as 
we are in a strong coupling regime and the calculation is non-perturbative. 

We cannot explicitly solve (3{a,B) = because we do not know the 
functional form of (3. This has to be done numerically |jT9| and will be 
reported elsewhere. Nevertheless one can make the following observation. 
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For a given amount of energy, if the dilaton has higher L it must take longer 
to reach the horizon. This is because higher L would mean its radial velocity 
is smaller and also the "centrifugal barrier" is higher. Thus thermalization 
should take longer. This in turn means that the Lyapunov exponents should 
decrease with increasing k"^. Thus for instance one can expect a relation of 
the form 

7(fc^) = 7o - cfc^; c > 0. 

Thus there should be a critical k"^ at which 7 = and there would be no 
thermalization. This corresponds to the impact parameter being larger than 
the black hole radius. But again as mentioned in the introduction, this is 
due to the neglect of the drag due to the background Hawking radiation. 
Presumably if this effect were included, 7, while decreasing with increasing 
fc^, would always be greater than zero. This seems to agree with numerical 
work on Yang- Mills done on the lattice |TT]. 



7 Conclusions 

We have discussed in this paper a possible resolution of a puzzle regarding 
black holes, which is simply stated: How does a pure state evolve into a state 
that appears to be at finite temperature? This question is difficult to answer 
directly in a theory of gravity. However using the AdS/CFT connection it 
can be posed as a problem in Yang-Mills. Namely if some energy is added in 
a coherent form say as an electric field excitation, at zero temperature, how 
does the theory "suddenly" acquire a finite temperature? That this must 
happen follows from the AdS/CFT connection. An answer to this question 
is that the theory "self-thermalizes" . This is known to be property of chaotic 
systems and Yang-Mills is classicaly known to be chaotic. This seems to 
answer the question qualitatively. 

Note that there is no implication here that quantum mechanics violates 
(or does not violate) unitarity. The statement is that the final state with the 
black hole appears to be at finite temperature. We are invoking the argument 
[0] that a pure state can have this property, in a chaotic system. 

A more quantitative check is to calculate the Lyapunov exponents in 
supergravity and compare with Yang-Mills calculation. An agreement would 
not only strengthen the AdS/CFT conjecture, it would also say that chaotic 
behaviour is there and can explain many otherwise mysterious facts about 
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black holes. We have set up the calculation and shown that the existence of 
a horizon implies that there are complex eigenvalues. The dependence of the 
Lyapunov exponent on temperature is as expected from Yang-Mills theory on 
dimensional grounds. The actual calculation involves numerical work. The 
infrared divergences that show up in perturbative QCD (for the non leading 
terms) do not show up here. We also have an equation for the dependence 
of the Lyapunaov exponent on /c^. But again solving it requires numerical 
work. The qualitative fact that the exponent should decrease with increasing 
k can be argued on physical grounds from the supergravity viewpoint. This 
is also in accord with |T0|, |1T[. 

Thus the role of chaos seems to be important. But more work remains to 
be done to quantitatively support this idea. 
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A Appendix: Illustrating Complex Bound- 
ary Conditions in Quantum Mechanics 

We illustrate how complex boundary conditions lead to complex eigenvalues 
in a simple one-dimensional quantum mechanical problem. Consider the 
potential shown in fig 2. 

There is an infinite barrier on the left and a finite area barrier on the 
right, which we represent by a (5-function of strength V for simplicity. We 
have a particle in the box that, for infinite V is confined inside. But for 
finite, large V will slowly leak out. The wave function is given by 

iP = A{k)e'^'' + B{k)e-'^'' : x < a. 

^ = C(fc)e''=" : X > a. (A.l) 

Note that outside the box we have assumed that the wave is purely right 
moving. This is a complex boundary condition (at infinity). This is analogous 
to the dilaton having only ingoing waves inside the horizon. 
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Fig 2. Particle leaks out of a box . Outside the box 
there is no left-moving wave and corresponds 
to the inside of a black-hole horizon. 



The other boundary conditions are 

ip{0) = 



(A.2) 



at X = and integrating Schroedinger's equation across the (5-function gives 
at X = a, in addition to continuity of ip, 



dip 
dx 



dip 
dx 



-Vi){a)2m. 



(A.3) 



We look for "stationary" solutions with a time dependence e with E 
possibly complex. Thus 

^ = E (A.4) 



2m 

Imposing ( [A .21 ) gives A = —B, and imposing continuity at x = a gives 



1 - e 



~2ika 



(A.5) 



Thus all constants A,B,C are determined upto an overall constant that 
can be fixed by normalization. However ( |A.^ ) still has to be satisfied and 
this gives an eigenvalue equation for k 



-2ika 



2ik- 



1-e 



-2ika 



-V. 



(A.6) 
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Let us set 

k = p + i n 

When V = oo the solution clearly is = and p = For large, finite 
V we can do a perturbation expansion in powers of ^ around this solution. 
We assume the wavenumber decreases a bit: 

P= ^ ;/i>0■ 
a(l + e) 

and e, /i are small and non-zero. 

Plugging into (|A.6|) we find to lowest order 



1 

^ ~ Va 

The fact that /i is negative means that the exponential is growing outside 
the box. In order that the wave function be normalizable, we have to assume 
that somewhere outside (but at a finite distance) there is a mechanism for ab- 
sorbing the particle. This would be a physical justification for the boundary 
condition. In the case of the black hole this is what we assume happens, once 
the particle enters the interior of the black hole and reaches the singularity. 

2 2 

Thus the real part of the energy is ^ ~^ and the imaginary part is the 

width r = ^ 

r gives the rate of escape of the particle from the box (=fiux) and can be 
seen to be equal to | C p x velocity, where velocity = ^. 

This example illustrates how energy becomes complex when complex 
boundary conditions are used to describe leakage of fiux. One can also 
check that the equation analogous to (|5.16|) for the Schroedinger equation 
is satisfied for this problem. 



References 

[1] J. D. Bekenstein, Phys. Rev D 7 (1973) 2333. 

[2] S. W. hawking. Comm. Math. Phys. 43 (1975) 199. 



17 



[3] M. Srednicki, Phys. Rev E 50 (1994) 888, |cond-mat/94U3U51| ; |cCT 
mat/94l"0046l 



J. Maldacena, Adv. Theor. Math. Phys. 2 (1998), |EeFth/9711200 



S.S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B428 
(1998) 105, |hep-th/9802109 . 



E. Witten, Adv. Theor. Math. Phys. 2 (1998), |hep^th/9802150 
E. Witten, |hep-th/9803T3T 



R. D. Pisarski, |hep-ph/9302242t ; Phs. Rev. Lett. 63, 1129 (1989). 

E. Braaten and R. D. Pisarski, Phys. Rev. Lett. 64,1338, (1990). 

T. S. Biro, S. G. Matinyan and B. Mueller, "Chaos and Gauge Field 
Theory", (1994),World Scientific, Singapore ; T. S. Biro, C. Gong and 
B. Mueller, Phys. Rev. D52 (1995) 1260, hep-ph/940932. 

B. Mueller and A. Trayanov, Phys. Rev. Lett. 68 (1992) 3387 ; T.S.Biro, 
M. Feurstein, H. Markum, |hep-lat/97lT002 . 

C. Csaki, H. Ooguri, Y. Oz and J. Terning, |hep-th/ 9806^021 . 

R. de Mello Koch, A. Jevicki, M. Mihailescu and J. P. Nunes, |hep- 



th/9806125 . 



M. Zyskin, |hep-th/9806T28 . 

T. Banks et al, |hep-th/9808016t V. Balasubramanaian, P. Kraus, A. 
Lawrence, S. Trivedi, Phys. Rev. D59 (1999)104021, |hep-th/98080T7 . 



D. Kabat and G. Lifshytz, JHEP 9812 (1998) 002, |hep-th/98062T4] , f^- 
th/9902073|; S. R. Das, JHEP 9902 (1999) 012, |hep-th/990lOTq, Eep 



th/9905037 



th/9901012 



G. Horowitz and N. Itzhaki, JHEP 9902 (1999) 010, lep 
D. Bak and S. J. Rey, |hep-th/990210ll . 



L. Dolan and R. Jackiw, Phys. Rev D9 (1974) 3320. 

H. A. Weldon, Phsy. Rev. D28 (1983) 2007. 

S. Bal, S. Kalyana Rama and B. Sathiapalan, work in progress. 



18 



